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The main thermodynamical properties of the first order phase transition of the relativistic mean- 
field (RMF) hadronic model were explored in the isobaric, the canonical and the grand canonical 
ensembles on the basis of the method of the thermodynamical potentials and their first derivatives. 
It was proved that the first order phase transition of the RMF model is the liquid-gas type one 
associated with the Gibbs free energy G. The thermodynamical potential G is the piecewise smooth 
function and its first order partial derivatives with respect to variables of state are the piecewise 
continuous functions. We have found that the energy in the caloric curve is discontinuous in the 
isobaric and the grand canonical ensembles at fixed values of the pressure and the chemical poten- 
tial, respectively, and it is continuous, i.e. it has no plateau, in the canonical and microcanonical 
ensembles at fixed values of baryon density, while the baryon density in the isotherms is discontin- 
uous in the isobaric and the canonical ensembles at fixed values of the temperature. The general 
criterion for the nuclear liquid-gas phase transition in the canonical ensemble was identified. 
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43 I- INTRODUCTION 

j The study of the phase transitions in the nuclear matter under extreme conditions is important for the theoret- 
ical analysis of heavy ion collision experiments at the intermediate and high energies [l|, Q and for understanding 
of many critical issues in astrophysics Q . It is generally accepted that the nuclear multifragmentation observed 
in intermediate-energy nuclear reactions is indicating the first order phase transition of the nuclear liquid-gas 
type [J] . This claim is effectively based on the concept of the plateau in the nuclear multifragmentation caloric 
curve (the dependence of the temperature T on the excitation energy E* of the system) which was theoretically 
predicted in [5[ and later was experimentally found in the 600 MeV/nucleon Au + Au collision Q. However, 
in certain experiments, as for example in 1 GeV/nucleon Au + C collision [tJ and 95 MeV/nucleon Ar + Ni 
collisions Q, the plateau in the experimental caloric curve is absent. Moreover, although there is a class of 
statistical multifragmentation models (SMM) which predicts a plateau (the discontinuity of energy) in the 
caloric curves in the canonical and microcanonical ensembles 0,111, 043 , there is yet another class of SMM 
which predicts a plateau in the caloric curves only in the isobaric ensemble and a very smooth transition in 
the canonical and microcanonical ensembles [15J420( . This particular property of the caloric curves for the first 
order phase transition in the canonical and microcanonical ensembles can be also seen in the models of other 
physical systems [2ll - [29l ]. 

The first order phase transition in almost all classification schemes has the Ehrcnfcst definition, which is one 
associated with a finite discontinuity in one or more of the first derivatives of the appropriate thermodynamic 
potential with respect to its variables of state [3(1 EJ . The concrete type of the first order phase transition 
is related to an appropriate potential f32j . For example, for the liquid-gas phase transition the Gibbs free 
energy G is relevant and there are discontinuities in its first derivatives, the entropy and the volume, across 
the temperature-pressure coexistence curve (3lj . Herewith, the energy in the caloric curve for the liquid- 
gas phase transition is discontinuous in the isobaric ensemble and it is continuous in the canonical and the 
microcanonical ensembles. Otherwise, for a magnetic system the free energy F is the appropriate potential with 
the discontinuity in its first derivative, the magnetization, across the temperature-magnetic field coexistence 
curve [32| . The energy in the caloric curve for the first order phase transition associated with the free energy F 
is discontinuous in the canonical and the microcanonical ensembles, as the first derivative of F, the entropy, is 
discontinuous [33| . Some other alternative methods which deal with the occurrence of phase transitions can be 
found in (H-E^ 

The nuclear liquid-gas phase transition, i.e. the first order phase transition associated with the Gibbs free 
energy G, usually occurs in the models of the interacting hadron gas such as the relativistic mean-field hadronic 
model [43T]49l |. the mean-field approximation with a Skyrme effective interaction [5(j, etc. In spite of the fact that 
these models have been intensively investigated during a long period of time, the liquid-gas phase transition in 
such models has not been studied completely. Therefore, the main purpose of this paper is to provide a holistic 
and concerted approach to tackle the problem of the nuclear liquid-gas phase transition on the basis of the 
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method of the thcrmodynamical potentials and their first derivatives in different statistical ensembles and to 
establish the general criteria for the liquid-gas phase transition which allow to distinguish it from the first order 
phase transition associated with the free energy F. To determine these criteria we use the well-known RMF 
hadronic model [43|, H|[ , which enables to find the consistent behavior of the caloric curve and the equation of 
state (the discontinuity of the density in the isotherms [5lL (HJ) for the nuclear liquid-gas phase transition in 
different statistical ensembles. 

The structure of the paper is as follows. In Section [Til wc briefly describe basic ingredients of the relativistic 
mean-field model and the methodology to evaluate the model parameters. The thermodynamic results for 
the liquid-gas phase transition and the caloric curve are discussed in Sections IIII1 The main conclusions are 
summarized in the final section. 



II. RELATIVISTIC MEAN-FIELD HADRONIC MODEL 



Let us consider the relativistic hadron many-body problem in the framework of the local quantum field theory. 
In the Lagrange formalism this theory is based on the action functional of the form S = J d 4 xC, where d 4 x is 
the integration measure in four-dimensional Minkowski space and L is the Lagrange density. We now proceed to 
use the simplest realization of the RMF approach [43[ . The effective Lagrangian of the system of the low-lying 
baryons interacting through the exchange of the scalar meson a and the intermediate abelian t/(l) gauge boson, 
i.e. vector meson to, is given by [iij 

C = l - ^Drf - (£>>)7 M V] " - -F^ + ^m 2 v A^ + ^(d^d^ - m 2 s 4> 2 ) + <? S #V, (1) 

where tp = ( , 2 ] is the isodoublet of the group SU(2), ip 1 and ?A 2 are the spinor fields for protons and neutrons 



with common mass m, 7'' is the Dirac gamma matrix, is the real isoscalar vector Proca field for the u meson 
with mass m v , <p is the real isoscalar scalar field for a meson with mass m s and g Vl 9s are the coupling constants, 
respectively. The antisymmetric field strength tensor for gauge field and the covariant derivative of the group 
U(l) are defined by 

F tu/ = d v ,A v -d v A ll , D /J , = d ll +tg v A^. (2) 

Note that throughout the paper we use the natural units system, K = c = feg = 1. 

This Lagrangian resembles the quantum electrodynamics (QED) Lagrangian with a massive abelian "photon" 
and an additional scalar interaction. This theory is renormalizable [4J, l53[ . The Lagrangian ([1]) is invariant 
with respect to the global gauge transformation of the abelian group U(l) and the non-abelian group SU(2). 
But, it is not invariant with respect to the local gauge transformations of these two groups. The mass term 
m^A^A^ in the Lagrangian ([1]) destroys the local abelian gauge U(l) invariancc. In order to be invariant with 
respect to the local gauge transformations of the non-abelian group SU(2), the Lagrangian ([TJ must contain an 
isotopic triplet of massless vector fields interacting with the spinor fields ip 1 and ip 2 , because the theory with 
non-abelian massive gauge fields is not renormalizable. If that is the case, to preserve the renormalizability of 
the Lagrangian, the mass terms for the non-abelian gauge quanta should be introduced into the theory via the 
so-called Higgs mechanism [53j | . However, in isospin symmetric matter the long-range part of the nuclear force 
which comes from one pion exchange and the force mediated by the exchange of the rho meson average to zero. 

In the mean-field approximation (44l |. the following conditions are imposed on the fields A^ and 0, 

= 0o, A ll = $A , (3) 

where the quantities 4>o and Aq are constants, independent of (i.e. d^A^ = and <9 M 0o = 0). In the second 
quantization, the energy operator, H = J d 3 xT 00 , the net baryon charge operator, B = J d 3 xip^f°ip, and the 
Euler-Lagrange equations for the classical fields Aq and 0o can be written as [44[ 

& = £ [ £ (+)°k a ifc + £ (-) b la h P<y] - \i m l A l - m 2 s4>l)V, (4) 

pa 

B = ^2 [dfcCifc - bi a b Pa .], (5) 

pa 

A ° ~ = °> ( 6 ) 

v pa 
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g a 1 ^ m t 



m V ^ E 

per 



(7) 



where £(±) = E* ± g v A is the zero component of the four momentum, P(±j = (£(±)jP)) °f the plane- wave 

solutions of the Dirac equation, E* = ^/p 2 + m* 2 is the zero component of the four momentum p^ = (E* , p) 
of the Dirac bispinor of amplitudes u± p ± a , m* = m — g s 4>Q, is the nucleon effective mass, p 2 = m* 2 and V is 
the volume of the system. We use the composite index notation, a = (s, J), where s and / are the nucleon spin 
and isospin components, respectively. The creation and annihilation operators for fermions, a~ a and ap a , and 

anti-fermions, bf a and &t , satisfy the equal-time anticommutation relations. Note that in Eqs. (0]) - (JT)) the 
vacuum terms were neglected. 

Let us consider a system of volume V, in contact with a heat and particle reservoir of temperature T and 
chemical potential p.. Then, the statistical operator and partition function are 
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y Z 



and 



Z = Tr e 



respectively. Plugging Eqs. (jU)-© into ©, we obtain \& 



-TlnZ 



e (+)->* e (-) +> * 
ln(l + e t ) + ln(l + e r ) 



pa 



H-)' 



V 



/As 



1 r— s m* 



per 



with the mean occupation numbers given by the Fcrmi-Dirac distribution functions 

1 , , 1 



{npa) 



{Ufa) 



e t 



1 



e t 



1 



and 



A (T) = ^( a vP 2 Ta s p 2 s )V, e {±) =E*±a v p, 



m =m - a s p s , 



(8) 

(9) 
(10) 

(11) 
(12) 
(13) 



where p = A g v /a v , p s = 4> g s /a s , a v = gl/m 2 , = C 2 /m 2 , a s = g 2 /m 2 s = C 2 /m 2 , with C v and C s being 
dimensionless parameters of the model [44j . 

Equations ((^|)— (fT5j) constitute a closed system of equations which allow us to determine all ensemble averages. 
Using Eqs. (|9j)- (|ll[) . we observe that the thermodynamic potential is independent of p and p s : 



~dp~) 



TVpp s 



vV(p - — ^[(ripv) - (n?*)}) = 0, 



pa 



dp s J 



a s V{p s -y^2 -^rl( n P<?) + (%*)]) = 0. 



In particular, the fundamental equations of thermodynamics are 

dQ = -SdT - pdV - (B)dfi, 
TdS = d(H) + pdV - fid(B) 



and the Eulcr theorem is 



TS=(H) + P V-p(B). 



(14) 
(15) 



(16) 
(17) 



(18) 



The mean energy (H), the mean net baryon charge (5), the pressure p, and the entropy 5* can be written as 

p(B) =Y,E*[(n fa ) + (n Pa )] + A (+) , (19) 



Vp. 
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(B) = =Y^M ( 20 ) 

>--{w)-&x&m+™ + ¥- (2i) 

s = -(§5) v =i(-n + (ff)-M<s». (22) 

where p = (i?) /V is the density of the net baryon charge. 

Let us introduce the notation e = {H)/V. The constants a s and a v are derived from the lowest-energy state 
of the system at T = and the conditions £b(po) = £bo &ud p = po, where = (e/p) — m is the binding energy 
per nucleon and po is the normal nuclear density. For the isospin symmetric nuclear matter Ebo = — 16 McV 
and po = 0-16 fm -3 , respectively (44|. We obtain (54[ 



1 ( Pfo 



ar 1 



e w + m-mjl + xZ), (23) 

3 ^2(^0) / PfoV 1 , 0/A 
h{x ) 31 2 {x ) ( Pfo\ Pfo 

— . 4 _ m — Ftn — m. 



4 pfo ^t + r " w " £b0 m x/ 1 + x ° (25) 



and 



& v J„(a:) = xx/l + x 2 (l + c v x 2 ) - ln(z + \/l + x 2 ), (26) 



where v = 1,2, 3, x- = Pfo/wq, = m*(po), Pfo = (Sir 2 Pa H) 1 ^ is the Fermi momentum, 7 is the spin-isospin 
degeneracy factor for nuclcons, 61 = 62 = 1, &3 = — 1 and ci = 2,C2 = 0,C3 = —2/3. The number Xq (or mg) 
is the solution of Eq. ([25]) . Substituting .To (or mj) into Eqs. ([23]) and (|24p . we obtain the numerical values for 
the parameters of the model at T = 0: C 2 ~ 250 and C 2 S ~ 330. 

III. FIRST ORDER PHASE TRANSITION FOR NUCLEAR MATTER 
A. Grand canonical ensemble (T, V,j£) 

Let us describe the phase transition of the RMF model in the grand canonical ensemble. This ensemble is 
characterized by the potential, f2(T, V, p) ([9]) and the physical observablcs correspond to its first and second 
order partial derivatives, Eqs. ([2T)|). (f2~Tj) and (|22p . The geometry of the hypersurface f2 determines the physical 
properties of the system. The thermodynamical potential (O is a homogeneous function of first degree of the 
extensive variable of state V, so in the follwing we shall work with the specific thermodynamic potential u> (the 
pressure p) , 

u(T,n) = Sl(T,V,n)/V=-p(T,n). (27) 
From (|27p we obtain the entropy and baryon densities, 

S fdp\ ._ , (B) (dp 



s{T ^^v = [w)^ ^ s ¥ = {t) T ' (28) 

which satisfy Eqs. ([15]). ([IT]), and (f!5]l. as 

tip = sg?T + pdp, Tds = de — pdp, Ts = e + p — pp. (29) 

The thermodynamic functions are multiple valued, as it can be seen in Figs. [T] [5] and [31 and therefore the 
stable states are obtained by the Maxwell construction, by excluding the loops-the areas A, B, and C in Fig. 
[1] The self-intersection points a on the curves p(p) and p(T) are found from the solutions of the equations 



Pi (p) = P3 (p) , Pi ( T ) = P3 ( T ) , 



(30) 
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FIG. 1: (Color online) The Maxwell construction for the specific thermodynamical potential p(T,p) (the pressure), its 
first order partial derivative, p = dp/dp (upper part) and for the pressure p as a function of the specific volume v (lower 
part) of the RMF model in the grand canonical ensemble at fixed temperature T (left panels) and at fixed chemical 
potential p (right panels). The curves depict the exact anomalous results (dotted lines) and the results with the Maxwell 
construction (solid lines). Symbols denote the phase transition point at T — 14.62 MeV and p = 914.95 MeV. The 
normal specific volume vo = 1/po = 6.25 fm 3 , where po is the normal nuclear density, po = 0.16 fm -3 . 



where p\ and p^ are the values of the specific thermodynamical potential p (the pressure) on the lines 1 and 3, 
respectively. These solutions are the numbers p*(T) and T*(p). In this way, the loops A on the curves p(p) 
and p(T) in Figs. [Ha) and[ljb), respectively, are replaced by the points a, whereas the beckbending curves 
connecting the points b and c on the functions p(p) and p(T) (Figs. [TJ; andQJl) are replaced by straight lines. 
This definition of the Maxwell construction is equivalent with the equal-area definition for p — v dependence 
[52| as represented in the the bottom panels of Fig. [TJ 
The critical point, (T c , p Cl p c ), is obtained from the equation 



mm 



dp 



(31) 



and for the typical numerical parameters of the RMF model, with a nucleon mass, m = 939 MeV, we find 
T c = 18.927 MeV, p c = 910.033 MeV and p c = 0.441p - 

Figure [2] represents the behavior of p, s, p, f, e and m*, as functions of the chemical potential p at fixed 
temperature T; f = F/V is the density of free energy. For T > T c , the thermodynamical variables p(p), p(p), 
s(p,), e(p), f{p) and m*(p) are continuous, one-valued, monotonic and differentiable functions of p (the line 
3,4 in all panels of Fig. [2|. For T < T c , the specific grand potential p(p), with the Maxwell construction, is 
a piecewise smooth function. The point of phase transition, denoted by p* , is a point of discontinuity of the 
first derivative of p (the solid lines 1,2 in Fig. [2]). For T < T c , with the Maxwell construction, p is continuous 
and single-valued for all p, but its first order partial derivatives with respect to variables of state, namely p(p) 
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FIG. 2: (Color online) Grand canonical ensemble (T,fi). The specific thermodynamical potential p (the pressure), the 
entropy density S/V, the baryon density p, the density of free energy F/V, the energy density {H)/V and the effective 
nucleon mass m* as functions of the chemical potential p at fixed temperature T for the RMF approach. The curves 
1, 2, 3, 4 are obtained at T = 7, 15 MeV, T = T c and T = 22 MeV, respectively. The solid lines 1, 2 are the results with 
the Maxwell construction. The symbol is the critical point and the dashed line is the phase diagram. 



and s(/ii), are discontinuous at /i = /i*. The jump of the entropy density at n = \x* is related to the latent 
heat-A = T5s. Moreover, p, p and s, are strictly increasing functions of fi. 

Because of the small values of p and the slow increase of p with p, we say that the region p < p* corresponds 
to the gas phase of the nuclear matter. 

In the region p > p* the baryon density is considerably larger and the presssure increases more rapidly with 
p. so we call this the liquid phase of the nuclear matter. 

The baryon density, the entropy density, the density of free energy and the energy density have positive 
jump discontinuities throughout the phase separation at the point of phase transition, p*, whereas the effective 
nucleon mass m* has a negative jump. 

In Fig. [3]we plot p, s, p, /, e and m* as functions of T, at fixed p. If p < p c or p > po (po is the chemical 
potential on the coexistence curve at T — 0-see Fig. 2]), the functions p(T), p{T), s(T), e(T), f(T) and m*(T) 
are continuous, single- valued, monotonic and differentiable for any T (the lines 1 and 2 in Fig. [3]). At fixed 
p c < p < po, the function p(T) with Maxwell construction is a piecewise smooth function with a point of 
discontinuity of the first order derivatives at T — T* (the line 3 in Fig. [3]). In this point p(T) is continuous, but 
has a sharp corner (cusp). The first order partial derivatives of p with respect to variables of state, namely p(T) 
and s(T), with Maxwell construction are single- valued piecewise continuous functions with jump discontinuities 
at T = T* . Therefore, according to the Ehrenfest clasification for phase transitions, the system undergoes a 
first order phase transition at T = T* , p being fixed in the interval p c < p < pq. 

The function p(T), p(T) and s(T), are strictly increasing functions of T. At T < T*, the baryon density is 
small and the system is in the gas phase, whereas at T > T* the baryon density is large and the system is in the 
liquid phase. The entropy density, the baryon density, the density of free energy and the energy density have 
positive jump discontinuities throughout the phase separation at the point of phase transition T*, whereas m* 
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FIG. 3: (Color online) Grand canonical ensemble (T,p). The specific thermodynamical potential p (the pressure), the 
entropy density S/V, the baryon density p, the density of free energy F/V, the energy density {H)/V and the effective 
nucleon mass m* as functions of the temperature T at fixed chemical potential p for the RMF approach. The curves 
1, 2, 3 are calculated for p = 905 MeV, p = p c = 910.033 MeV and p — 914.95 MeV, respectively. The continuous curve 
3 is the results with the Maxwell construction. Symbol depicts the critical point. The dashed line is the phase diagram. 

has a negative jump discontinuity. 

The phase diagrams for the symmetric nuclear matter in the grand canonical ensemble (T, p) in the RMF 
approach are depicted in Fig. |4j The phase diagram T — p for the first order (liquid-gas) phase transition is 
represented by the continuous coexistence curve which crosses the p axes at the point (T = 0, p, = po) and it is 
finished at the critical point (T c , p c ), where po = m + Sbo = 923 MeV. Along this curve, p(T, p) and s(T,p), are 
discontinuous. The phase diagrams T — p and p — p also are represented by the continuous coexistence curves 
because the function p(T, p) is continuous in the points of phase transition. However, in the phase diagrams 
T — p, T — s, n — p, and /i-swc have coexistence regions instead of coexistence lines, because the variables 
p and s are undefined in the points of phase transition (T* , p* ) . Note that the mixed phase is defined by the 
coexistence lines on the phase diagrams T — p, T — p and the coexistence areas on the phase diagrams T — p, 
T-s, 

In conclusion, the phase transition that appears in the nuclear matter in the RMF model in the grand canonical 
ensemble is accompanied by discontinuities in the first order derivatives of the thermodynamic potential CI and 
therefore is a first order phase transition in the Ehrenfest classification [30| . 

B. Isobaric ensemble (T,p,B) 

Let us investigate the phase transition of the RMF model in the isobaric ensemble. The thermodynamical 
potential of the isobaric ensemble is the Gibbs free energy G = Cl+pV + pB and its first order partial derivatives 
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FIG. 4: (Color online) The phase diagrams T — p, T — p, T — s, T — p (top panels) and p, — p, p — s, p — p (bottom 
panels) for the first order phase transition of the RMF approach. Roman numerals denote the homogeneous phases, 
the gas phase (7) and the liquid phase (//), the shaded areas correspond to the mixed phase, whereas the lines are the 
coexistence curves. Symbol depicts the critical point. 



with respect to the variables of state, (T,p, -B), are the entropy, the volume and the chemical potential, 

f9G\ (dG\ fdG\ 

S = -KSr) rB ' "-(*)„■ "-{sbK- (32) 

The potential G is proportional to -B, and therefore the specific thermodynamic potential is 

g(T,p) = G(T,p,B)/B = p(T,p). (33) 

The first order partial derivatives of the specific thcrmodynamical potential fi(T,p) with respect to the variables 
of state (T, p) can be written as 

«*■'>--(£),■ »<^KtX- <-» 

where s — S/B is the entropy per nucleon and v = V/B is the specific volume. They satisfy the relations 

dp, = —sdT + vdp, Ts = e + pv — p, (35) 

where e = (H) / B is the energy per nucleon. 

In Fig. [5] we plot the functions fj,, s, v, /, e and m* /m as functions of p, at fixed T for the RMF approach; 
/ = F/B is the free energy per nucleon. For T > T c , the functions p(p), v(p), s(p), e(p), f(p) and m*(p) are 
continuous, one-valued, monotonic and diffcrentiable for any p (the lines 3, 4 in Fig. |SJ). For T < T c , p(p) is still 
a continuous function, but its first order derivatives are discontinuous at the point of phase transition, p* (the 
solid lines 1 and 2 in Fig. [5j. The first order partial derivatives of the thermodynamic potential with respect to 
variables of state, i.e. the functions v(p) and s(p), with the Maxwell construction, have jump discontinuities at 
p — p* . Therefore in the Ehrenfest classification this is also a first order, liquid-gas, phase transition. 
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FIG. 5: (Color online) Isobaric ensemble (T,p). The Gibbs free energy per micleon p (the chemical potential), the 
entropy per nucleon s, the specific volume v, the free energy per nucleon /, the mean energy per nucleon, e and the 
effective nucleon mass m* as functions of the pressure p at fixed temperature T for the RMF model. The curves 1, 2, 3, 
and 4 are obtained for T = 7, 15 MeV, T — T c and T = 22 MeV, respectively. The solid lines 1 and 2 correspond to the 
Maxwell construction. The symbol is the critical point and the dashed line is the phase diagram. 



The jump of the entropy per nucleon at the phase transition pressure, p*, is related to the latent heat. The 
functions [i{p) and f(p) are increasing with p, whereas v(p), s(p), e(p) and m*(jp) are decreasing. In the gas 
phase, p < p* , the specific volume, v(p), the entropy per nucleon, s(p), and the energy per nucleon, e(j>), have 
large values and the reduced effective nucleon mass is close to one (m*(p)/m ~ 1). In the liquid phase (p > p*), 
the functions v{jp), s(p), e(p) and m*(p)/m take smaller values. The functions v(p), s(p), e(p) and m*{p) have 
negative jumps at the phase transition pressure, p* , whereas the free energy per nucleon f(p) has the positive 
jump discontinuity. 

Figure [S] presents the behavior of p, s, v, /, e and m*/m, as functions of temperature T at fixed p for the 
RMF approach. At the pressures higher than p c , the functions n(T), v(T), s(T), s(T), f(T) and m*(T) are 
continuous, single- valued, monotonic and diffcrcntiable functions on all of T (the lines 3 and 4 in Fig. [5]). 

At pressures below p c , the specific Gibbs potential, p(T), with the Maxwell construction, is a piecewise 
smooth function with a singularity at the phase transition temperature, T* (the solid lines 1 and 2 in Fig. [6]). 
The first order partial derivatives of /i with respect to variables of state, i.e. the functions v(T) and s(T), are 
single- valued piecewise continuous and smooth function for all T > 0, except in T = T* , where they have jump 
discontinuities. Also the functions /, i and m* have jump discontinuities at the point of the phase transition. 

The discontinuities of the quantities s(T), v(T), e(T) and m*(T) are positive, whereas the discontinuity of 
f(T) is negative. Therefore in this ensemble also the phase transition is of the first order, according to the 
Ehrenfest clasification. 

The functions p(T) and f(T) are monotonically decreasing, whereas the functions v(T), s(T), e(T) and m*(T) 
are monotonically increasing. In the gas phase, which correspond to a temperature T > T* , the specific volume 
v(T), the entropy per nucleon s(T) and the energy per nucleon e(T) have large values, whereas the reduced 
effective nucleon mass is close to one (m*(T)/m ~ 1). In the liquid phase, T < T* , the functions v(T), s(T), 
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FIG. 6: (Color online) Isobaric ensemble (T,p). The Gibbs free energy per nucleon, /j, (the chemical potential), the 
entropy per nucleon, s, the specific volume, v, the free energy per nucleon, /, the mean energy per nucleon, e, and the 
effective nucleon mass, m* , as functions of T at fixed p for the RMF approach. The curves 1, 2, 3, and 4 are obtained 
at the pressure p = 0.05 MeVfm~ 3 ,0.2 MeV fm -3 , p = p c = 0.484 MeV fm~ 3 and p = f.O MeV fm -3 , respectively. The 
solid lines f and 2 are the results of the Maxwell construction. The symbol is the critical point at p c = 0.484 MeV fm~ 3 . 
The dashed line is the phase diagram. 



e(T) and m*{T) take smaller values. 

Summarizing, we showed that in the isobaric ensemble the thermodynamic potential, fj,(T,p) (the chemical 
potential), is a continuous function, with piecewise continuous derivatives. The first order derivatives with 
respect to T and p, which are the entropy per nucleon s and the specific volume v, respectively, have jump 
discontinuities at the phase transition. The phase diagram in the coordinates T— p is represented by a continuous 
coexistence curve (Fig. @Jl) , whereas in the coordinates T — v and T — s is represented by the coexistence regions 
because the variables v and s are undefined in the points of phase transition at fixed values of T and p. 

Being accompanied by a jump discontinuity of the first order derivatives of the thermodynamic potential G, 
the phase transition of the RMF model in the isobaric ensemble is a first order phase transition in the Ehrenfest 
classification [30| and, therefore, it is a liquid-gas type [3lL l32|. 



C. Canonical ensemble (T, V, B) 

Let us investigate the phase transition of the RMF model in the canonical ensemble. The thermodynamic 
potential of the canonical ensemble is the Helmholtz free energy, F(T, V, B) = + fJ.B. The first order partial 
derivatives of F with respect to the variables of state, (T, V, B), are 

s= -(^L' p= -(^) TB > ^SL (36) 
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FIG. 7: (Color online) Canonical ensemble (T,v). The free energy per nucleon /, the entropy per nucleon s, the pressure 
p, the chemical potential /i, the mean energy per nucleon e and the effective nucleon mass m* as functions of the specific 
volume v at fixed temperature T for the RMF approach. The curves 1, 2,3,4 are obtained at T = 7, 15 MeV, T — T c 
and T — 22 MeV, respectively. The continuous curves 1,2 are the results with the Maxwell construction. The symbol 
depicts the critical point and the dashed line is the phase diagram. 



Since F is a homogeneous function in the variables V and B, in the following we shall work with the specific 
thermo dynamical potential (the free energy per nucleon) 

f(T,v) = F(T,V,B)/B. (37) 

The first order partial derivatives of the specific thermodynamical potential f(T, v) with respect to the variables 
of state (T, v) can be written as 

»P» = -(|) ; «?•<■■> — (%) t - < 38 > 

Due to the homogeneity property of the thermodynamical potential (|37|) all specific functions of variables of 
state of the canonical ensemble are intensive quantities. These quantities satisfy the differential equation for / 
and the Euler theorem 

df = —sdT — pdv, Ts — e + pv — /i. (39) 

To obtain the functions of the canonical ensemble we change the variables of state of the grand canonical 
ensemble (T, fi) to the variables (T, v) together with the Legendre transformation for the thermodynamical 
potential, f = pv — fx. 

Figure [7] presents the behavior of the specific thermodynamical potential /, the entropy per nucleon s, the 
pressure p, the chemical potential fi, the energy per nucleon e and the effective nucleon mass m* as functions of 
the specific volume v at fixed temperature T for the RMF approach in the canonical ensemble. At temperature 
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FIG. 8: (Color online) Canonical ensemble (T, v). The free energy per nucleon /, the entropy per nucleon s, the pressure p, 
the chemical potential p, the mean energy per nucleon e and the effective nucleon mass m* as functions of the temperature 
T for the fixed specific volume v (the baryon density p = 1/v) for the RMF approach. The curves 1,2,3 depict the 
results at the specific volume v/vo = 10 (p/po = 0.1), v = v c (p = p c ) and v/vo — 1.064 (p/po — 0.94), respectively. The 
solid lines 1, 2,3 are the results with the Maxwell construction. Symbol is the critical point at v c — 2.268uo- The dashed 
line is the phase diagram. 



T > T c , the exact thermodynamical quantities f{v), s(v), p(v), /x(i>), and m*(v) are continuous, one-valued, 
monotonic and differentiable functions on all of v (the lines 3,4 in Fig. [7]). At the temperature T < T c the 
specific Helmholtz potential f(v) with the Maxwell construction, its first derivatives, the entropy per nucleon 
s(v) and the pressure p(v), and the quantities n(v), e(v) and m*{v) are characterized by the linear changes in 
a closed interval vu < v < vi, where vj (vu) are the specific volume v in the gas (liquid) phase across the 
phase boundary at fixed temperature T (see the solid lines 1,2 in Fig. [7|). The chemical potential p,(v) and 
the pressure p(v) (isotherms) have the plateau as functions of v at fixed temperature T because they take the 
constant values in the points of phase transition. The function f(v) has the linear decrease, but the functions 
s(u), e(v) and m*(u) have linear growth in the region of mixed phase, vu < v < V[. In the liquid phase, at 
v < vji and T < T c , the entropy per nucleon s(v) and the effective nucleon mass m* decrease with decreasing 
of v, but the functions p(v), fi(v) and f(v) increase. The energy per nucleon e achieves his minimum at v = vq. 
In the gas phase, at v > vi and T < T c , the entropy per nucleon s(v), the energy per nucleon e and the effective 
nucleon mass m* increase with increasing v, but the functions p(v), p:(v) and f(v) decrease. 

Figure [5] presents the behavior of the specific thermodynamical potential /, the entropy per nucleon s, the 
pressure p, the chemical potential /x, the energy per nucleon e and the effective nucleon mass m* as functions 
of the temperature T at fixed specific volume v for the RMF approach in the canonical ensemble. The Maxwell 
construction at fixed values of v for any function A(T, v) is defined by the phenomenological equation 



A = Ai + a(Aji - Aj) 



P~ Pi 
Pi i - Pi 



(40) 



where Aj {An) and pi (pn) are the thermodynamical quantity A(T,p,) and the baryon density p(T,p:) in the 
gas (liquid) phase across the phase boundary in the grand canonical ensemble. Then, the function A(T, v) of 
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the canonical ensemble is obtained from the Eq. (|40|) by changing the variables of state (T, p) into the variables 
of state (T,v). At the fixed specific volume from the interval v < v < oo, the specific Helmholtz potential, 
f(T), with the Maxwell construction, its first derivatives, the entropy per nucleoli s(T) and the pressure p(T), 
and the quantities p{T), e(T) and m*(T) begin to deviate from their exact values at the temperatures T < T s , 
where T s is temperature of the crossing point on the coexistence curve. The first derivatives, s(T) andp(T), and 
the quantities /i(T) and e(T) are single-valued continuous broken-line functions and they have sharp corners 
(cusps) at the temperature T — T s . This point determines the transition of the system from the mixed phase 
at T < T s to the gas (liquid) phase at T > T s . Therefore, at fixed specific volume v the system can be only in 
two phase: mixed and gas (liquid). The functions p(T) and f(T) are strictly decreasing functions on a set /, 
but the functions s(T), p(T) and e(T) are strictly increasing functions on a set /. The functions p{T) and p(T) 
in the mixed phase at T < T s coincide with their coexistence curves because they take the constant values in 
the points of phase transition. 

Summarizing, we have obtained that in the canonical ensemble the first order phase transition of the nuclear 
liquid-gas type associated with the Gibbs free energy is defined by the linear changes of the specific thermo- 
dynamical potential / and its first order partial derivatives, the entropy per nucleon s and the pressure p, as 
functions of the specific volume v at fixed temperature T. Otherwise, it is also defined by the curvilinearly 
varying of the specific thcrmodynamical potential, the entropy per nucleon and the pressure, as functions of the 
temperature T at fixed specific volume v and by the sharp corners of the first derivatives s and p in the points 
where they cross the coexistence curves. The phase diagram T — v is depicted by the coexistence aria instead 
of the coexistence line. 



D. Caloric curve and the equation of state 

Let us investigate the caloric curve, i.e. the dependence of the temperature T on the energy E of the system, 
for the phase transition of the RMF model in the grand canonical, isobaric, canonical and microcanonical 
ensembles. The thcrmodynamical potential of the microcanonical ensemble is the entropy S(E, V, B). The first 
order partial derivatives of the thcrmodynamical potential S with respect to the variables of state (E, V, B) are 
the temperature, the pressure and the chemical potential 



1 (8S\ p (dS\ n fdS 



T \dE) VB T \8VJ EB T \dB /EV 



(41) 



The Legendre transformation for the thermodynamical potential is the equation TS = E — Q — pB. 

The thermodynamical potential of the microcanonical ensemble, the entropy S, is a homogeneous function of 
first degree of the extensive variables of state V, B 

S(E,V,B) = Bs(e,v), (42) 

where s(e,v) is the specific thermodynamical potential, v = V/B and e = E/B. Therefore, the variable B 
is excluded and the microcanonical ensemble is described by the specific thermodynamical potential s(i, v) 
as a function of the intensive variables of state (e,v). Then, the first order partial derivatives of the specific 
thermodynamical potential s(e, v) with respect to the variables of state (e, v) can be written as 

ne»=(g)~\ »(«»=<f;x. <«> 

Due to the homogeneity property of the thermodynamical potential (|42|l all specific functions of variables of 
state of the microcanonical ensemble are intensive quantities. These quantities satisfy the differential equation 
for s and the Euler theorem 

Tds = de + pdv, Ts = e + pv - fi, (44) 

To obtain the functions of the microcanonical ensemble we change the variables of state of the grand canonical 
ensemble (T, fx) to the variables (e, v) together with the Legendre transformation for the thcrmodynamical 
potential, s = (e + pv — /i) jT . 

Figure [9] presents the caloric curve in the isobaric, canonical, microcanonical and grand canonical ensembles, 
and the equation of state or the isotherms, i.e. the dependence p — p at fixed temperature T, in the canonical 
and isobaric ensembles for the nuclear liquid-gas phase transition of the RMF model. The excitation energy 
per nucleon is defined by the equation, E* = e(T) — e(0), where e(0) = 923 MeV is the energy per nucleon in 
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FIG. 9: (Color online) The temperature T as function of the excitation energy per nucleon E* (the caloric curve) for 
the RMF model in the different statistical ensembles at three points of phase transition (lines 1,2,3): (a) the isobaric 
ensemble at fixed pressure p, (b) the canonical and microcanonical ensembles at fixed baryon density p and (c) the grand 
canonical ensemble at fixed chemical potential p. (d) The isotherms in the canonical and isobaric ensembles at fixed 
temperature T (solid lines) in the same points of phase transition. Symbol is the critical point and the dotted line is the 
phase diagram. For details see the text and the Table U 



TABLE I: The parameters of the three points of phase transition for the RMF model 
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the ground state. The lines 1,2,3 in Fig. [9] were calculated for three points of phase transition. In Table [H 
the parameters of these points are summarized. In the point of phase transition at temperature T = T* the 
excitation energy per nucleon E* in the caloric curve has a jump discontinuity in the isobaric and grand canonical 
ensembles at fixed pressure p and at fixed chemical potential p, respectively. Otherwise, in the canonical and 
microcanonical ensembles at fixed specific volume v or baryon density p the excitation energy per nucleon E* in 
the caloric curve does not have a jump discontinuity. The curve T(E*) is a single- valued continuous broken-line 
function and has a sharp corner at the temperature T = T s , where it crosses the coexistence curve. In the 
mixed phase at T < T s the caloric curve T(E*) in the canonical and microcanonical ensembles is a continuous 
increasing function. 

As T increases, the system in the isobaric ensemble at fixed pressure p exhibits the phase transition from the 
liquid phase to the gas phase through an intermediate mixed phase created in the point of phase transition at, 
temperature T = T*. In the grand canonical ensemble at fixed chemical potential p the system exhibits the 
phase transition from the gas phase to the liquid phase also through the mixed phase. But, in the canonical and 
microcanonical ensembles at fixed baryon density p, as T increases, only the transition from the mixed phase 
to the gas (liquid) phase takes place. Moreover, at small temperatures and < p/po < 1 the nuclear matter in 
the canonical and microcanonical ensembles is situated in the mixed phase. In particular, the liquid-gas phase 
transition is characterized by the jump discontinuity of the baryon density p or specific volume v in the points 
of phase transition, i.e. it is defined by the plateau in the isotherms in the isobaric and canonical ensembles 
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at fixed temperature T. See Fig. [9] The uncertainty of p in the points of liquid-gas phase transition excludes 
the discontinuity of the energy in the canonical and microcanonical ensembles at fixed values of p, because for 
the one fixed value of the baryon density at constant temperature we enable to choose only one value of the 
energy. Therefore, the caloric curve for the liquid-gas phase transition in these two ensembles do not contain 
the plateau. 

Summarizing, we have found that for the nuclear liquid-gas phase transition of the RMF model the energy 
in the caloric curve is discontinuous in the isobaric and the grand canonical ensembles at fixed values of the 
pressure and the chemical potential, respectively, and it is continuous, i.e. it has no plateau, in the canonical 
and the microcanonical ensembles at fixed values of the baryon density. However, the baryon density in the 
equation of state (the isotherm) is discontinuous in the isobaric and the canonical ensembles at fixed values 
of the temperature. For one of the variants of the SMM the similar results were obtained in [17j . Thus, the 
general criterion for the nuclear liquid-gas phase transition, i.e. the first order phase transition associated with 
the Gibbs free energy G, in the canonical ensemble requires that the baryon density in the isotherms should 
be discontinuous at fixed values of the temperature and the energy in the caloric curves should be continuous, 
i.e. it should not have plateau, at fixed values of the baryon density (the baryon charge and the volume). Note 
that if for the certain physical system the energy in the caloric curve in the canonical and the microcanonical 
ensembles at constant values of the variables of state of the canonical ensemble is discontinuous, then its phase 
transition is a first order phase transition associated with the free energy F, because the discontinuity of the 
energy E, when F is continuous, is related to the discontinuity of the entropy, E = F + TS, which is the first 
derivative of the potential F [33||. We should also mention that all results of this work for the RMF model 
concerning the properties of the first order phase transition of the liquid-gas type associated with the Gibbs 
free energy G are in agreement with the principles of the general theory of phase transitions [13, H3, [32| . 

IV. CONCLUSIONS 

The first order phase transition for the RMF model was investigated on the basis of the method of the ther- 
modynamical potentials and their first derivatives in different statistical ensembles. The main thermodynamical 
properties of this phase transition were found by using the Maxwell construction in the framework of the grand 
canonical, canonical and isobaric ensembles. It was established that the first order phase transition of the RMF 
model is the phase transition of the nuclear liquid-gas type which can be associated with the Gibbs free energy 
G and its properties totally satisfy the requirements of the general theory of phase transitions. Indeed, we have 
found that in the isobaric ensemble the Gibbs free energy per nucleon (the chemical potential) for the RMF 
model is the piecewise smooth function and its first order partial derivatives with respect to variables of state 
(T,p), i.e., the entropy per nucleon and the specific volume, are the piecewise continuous functions. In the 
points of phase transition the chemical potential is a continuous function which has a cusp both as a function 
of p at fixed T and as a function of T at fixed p and the first order partial derivatives of the chemical potential 
with respect to variables of state, the entropy per nucleon and the specific volume, have jump discontinuities. 

Also, we have revealed that for the RMF model in the grand canonical ensemble the specific grand potential 
(the pressure) is the piecewise smooth function and its first order partial derivatives with respect to variables 
of state {T,p), i.e. the entropy density and the baryon density, are the piecewise continuous functions. In 
the points of phase transition the pressure is a continuous function which has a sharp corner (cusp) both as a 
function of p, at fixed T and as a function of T at fixed p, and the first derivatives, the entropy density and 
the baryon density, have jump discontinuities. However, in the canonical ensemble the definition of the nuclear 
liquid-gas phase transition appears unable to meet the criteria established for the grand canonical and isobaric 
ensembles. The first order derivatives of the Hclmholtz free energy per nucleon with respect to variables of state, 
the entropy per nucleon and the pressure, as functions of specific volume (baryon density) at fixed temperature 
vary linearly in the region of phase transition, but as functions of the temperature at fixed specific volume vary 
smoothly and have sharp corners in the points where they cross the coexistence curves. The first derivatives 
of the free energy have no jump discontinuities in the canonical ensemble. Thus, the nuclear liquid-gas phase 
transition in the RMF model is characterized by the jump discontinuities of the baryon density (the specific 
volume) and the entropy density (the entropy per nucleon) at constant values of temperature, chemical potential 
and pressure. This implies that for the nuclear liquid-gas phase transition of this model the phase diagrams 
T — p and T — p are represented by coexistence lines, however, the phase diagrams T — s (T — s) and T — p 
(T — v) are depicted by the coexistence areas. 

The caloric curve and the equation of state for the RMF model were calculated. It was established that 
the energy in the caloric curve for the nuclear liquid-gas phase transition of the RMF model is discontinuous 
in the isobaric and the grand canonical ensembles at fixed values of the pressure and the chemical potential, 
respectively, and it is continuous in the canonical and microcanonical ensembles at fixed values of baryon density 
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or the specific volume. However, the baryon density in the isotherms is discontinuous in the isobaric and the 
canonical ensembles at fixed values of the temperature. The general criterion for the nuclear liquid-gas phase 
transition in the canonical ensemble was identified. It states that the baryon density in the isotherms should 
be discontinuous at fixed values of the temperature and the energy in the caloric curves should be continuous 
at fixed values of the baryon density. It should be also mentioned that the obtained results for the first order 
phase transition of the RMF model are in total agreement with the principles of the general theory of phase 
transitions. 
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